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® Method and apparatus for communication channel identification and signal restoration by sending 


a test signal used for channel impulse response calculation. 


© A test signal is transmitted over a communication channel and the signal is received by a receiver. The 



temporal evolution of the test signal is represented in a matrix and the inverse of the test signal evolution matrix 
is stored in the receiver. Samples of the received signal and the test signal evolution matrix inverse are used to 
compute the channel impulse response, and the values of the channel impulse response set signal restoration 
filter coefficient values. The received signal is then restored by filtering. 
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BACKGROUND OF THE INVENTION 

The present invention relates to a method and apparatus for identifying the characteristics of a 
communication channel, and to a method and apparatus for restoring signals altered by the communication 
5 channel by using the channel characteristics information. 

Communication engineering continually must deal with the problem of restoring a signal which has been 
altered by the communication channel over which the signal was transmitted. Signal restoration often can 
be achieved if the communication channel is fully characterized, at least as to those parameters which 
contribute to the signal alteration. Thus, a frequently essential component of the signal restoration problem 
w is that of identifying the characteristics of the communication channel. 

A straightforward approach to the channel identification problem is to transmit a known signal over the 
channel, and to receive the transmitted signal after it has passed through the channel. The originally 
transmitted signal is compared with the received signal, and a model of the channel characteristics is 
developed based on the comparison. 
J5 If channel identification is to be carried out in real; time, or nearly in near real time, the computational 
burden can become enormous. In many practical applications hardware and cost constraints require 
considerable compromise as to the channel identification scheme that will be implemented, and may 
preclude anything approaching real-time channel identification. 

One area where hardware and cost constraints on signal restoration impose substantial limitations is 
20 that of echo cancellation, and in particular, ghost cancellation in video signals. A good survey of the ghost 
cancellation problem can be found in W. Ciciora et al, A Tutorial On Ghost Cancelling In Television 
Systems, IEEE Trans. On Consumer Elec., Vol. CE-25, No. t, pp 9-44 (Feb. 1979). 

It is an object of the present invention to provide a method and apparatus for restoring a received signal 
that has undergone alteration after passing through a communication channel. 
25 It is another object of the invention to provide a method and apparatus for identifying the characteristics 
of a communication channel. 

It is another object of the invention to provide a method and apparatus for restoring a received signal to 
eliminate echoes or video ghosts, as well as other signal distortions which may be present. 

It is another object of the invention to provide a method and apparatus for restoring a received signal 
30 which is simple and inexpensive, and which can be implemented with a minimum of computing hardware. 

Another object of the invention is to provide a method and apparatus for restoring a received signal 
which can operate in substantially real time in a television receiver. 

SUMMARY OF THE INVENTION 

35 

In the method according to the invention, no assumptions are made about the communication channel 
characteristics other than that the signal and channel interaction is described by linear system theory. 
Consequently, the channel is characterized completely by its impulse response. 

A test signal is transmitted over the channel and the received test signal will be distorted or otherwise 
40 altered by the channel. The test signal and the received signal are related by the channel impulse 
response. 

According to the invention, samples of the test signal are represented by an ensemble of sample 
sequences ordered in a matrix, called the test signal evolution matrix. The test signal evolution matrix 
changes from row to row in the manner that the test signal evolves in time. Samples of the received signal 
45 and the channel impulse response can be ordered in vectors. The channel impulse response sample vector 
can then be directly obtained through the multiplication of the inverse test signal evolution matrix and the 
received signal sample vector. 

The sequence of channel impulse response function samples thus obtained serves to provide the 
correct coefficient values to a filter. The received signal is then fed into the filter, and the received signal is 
50 restored by filtering. 

This method is general in nature, and is not derived from any particular model of signal distortion. In 
particular, it does not posit any particular model for echo generation such as the simple delay and 
attenuation model which is frequently assumed. 

Because the test signal is known, the test signal evolution matrix and its inverse are known. At least 
55 some of the values of the inverse matrix elements are stored in the receiver, and the multiplication of the 
inverse test signal evolution matrix and the received signal vector is carried out at the receiver using 
relatively simple and inexpensive hardware. 

When the method is applied to television reception, the test signal is contained in the vertical interval of 
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the transmitted video signal. The received test signal is processed according to the invention and the 
channel characteristics thus obtained are used to set the coefficient values of a video filter circuit. The filter 
coefficient values are set during the vertical interval so that the channel information is continually updated 
during the course of signal reception. 

5 

BRIEF DESCRIPTION OF THE DRAWING 

Fig. 1 is a schematic diagram of a circuit for computing values of a communication channel impulse 
response; 

10 Fig. 2 is a schematic diagram of another circuit for computing values of a communication channel 
impulse response; 

Fig. 3 is a graph of a test signal used in carrying out the invention; 

Figs. 4A and 4B illustrate how the test signal evolution matrix is written in cyclic Toeplit2 form; and 
Fig. 5 is a block diagram of a filtering arrangement for use in carrying out the invention. 

15 

DETAILED DESCRIPTION OF THE INVENTION 
Channel Identification 

20 The process of channel identification is based on a linear model for the interaction of the channel and 
the transmitted signal so that the received signal is described by the convolution 



n 

2s K(n)=£T(i)C(Ji-i+l) f 1 * 

where {T(i)} is a sequence of samples of a transmitted signal T, {C(i)} is a sequence of samples of the 
30 channel impulse response C, and {R(i)} is a sequence of samples of a received signal R which is the 
transmitted signal after it has passed through the channel. Throughout this description the indices i, j, k, L, 
M, N, etc. are integers with lower case indicating a running index and upper case indicating a particular 
value. 

Both of the sequences {T(i)} and {R(i)} are accessible to the television receiver if the signal T is known 
35 beforehand, e.g. if it is a predetermined test signal. The problem of channel identification then becomes one 
of calculating the sequence {C(i)} of channel impulse response samples from the sequences {T(i)} and {R- 

The invention makes the calculation of the sequence {C(i)} by expanding the sequence {T(i)} into a 
test signal evolution matrix [T(m,i)] (or simply [T]) which embodies information about the temporal evolution 
40 of {T(i)}. 
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75 



Each row "m" of the matrix is a sub-sequence of {T(i)} starting with the initial sample T(1) and having 
"m" members of the sequence, with each subsequent row having the next sample of {T(i)} appended to the 
sub-sequence of the previous row. If the sequences {C(i)} and {R(i)} are represented by respective vectors 
C(i)] and R(i)],then the convolution equation (1) defining the received signal can be written as 

25 

R(m)] = [T(m j)] C(i)] (2) 

This equation can be solved directly for the channel impulse response samples by multiplying both 
sides by the inverse of [T(m,i)] which yields 

30 

C(p)] = [T(p,q)]- 1 R(q)] (3) 

The matrix [T(m,i)] is comprised of samples of the test signal T and thus it is known beforehand and 
wholly within the control of the system designer. Consequently, so is the inverse matrix [T(p,q)]~ 1 . The 
35 elements of [T]~ 1 can be stored in the television receiver and the matrix multiplication carried out by 
successive convolutions as the sequence {R(i)} is being received, as follows: 

C(1) = T-'O, 1) R0) + T" 1 (1,2) R(2) + T" 1 <1,3) R(3) ... 
C(2) = T-^2, 1) R(1) + T~ 1 (2, 2) R(2) + ^(2,3) R(3) ... 

40 

Q 

C(p)=^T' l (p f q)R(q) 

45 

where T~ 1 (p, q) = is the (p, q)th element of [T(p, q)]~\ 

Computation of the matrix equation 3 can be simplified by adding additional elements to appropriate 
rows of [T] to create a new matrix [T] which is a cyclic Toeplitz matrix whose elements on each diagonal 
are identical and whose rows and columns are cyclic. A property of the cyclic Toeplitz form is that the 

so inverse matrix [T]~ 1 is also a Toeplitz form. Thus, a single row of the matrix can generate every other 
element of the matrix simply by shifting the second (column) index of values of the matrix elements of a 
particular row. The symmetrical property of the Toeplitz form allows elements from a single row of the 
inverse matrix to be used to carry out all of the convolution summations represented by equation (3). The 
convolutions which generate the sequence {C(i)} can be carried out with the circuit shown in Fig. 1 . 

55 In Fig. 1, a register 1 stores the sequence {T~ 1 (M, j) where M denotes a particular row of the matrix [T]- 
" 1 which contains N elements. The received signal sample sequence {R(i)} is stored in the shift register 2. 
The content of the respective cells of register 1 and shift register 2 are multiplied by the plurality of 
multipliers 3 and the products are added by adder 4. The output of the adder 4 is a value C(i) of the 
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channel impulse response sample sequence. The entire sequence {C(i)} can be generated by storing the 
first received signal sample R(1) in the shift register 2 and calculating C(1); shifting R(1) by one cell in the 
shift register 2, storing R(2) and computing C(2); and so forth. 

Another circuit for carrying out the calculations is illustrated in Fig. 2. The read-only memory (ROM) 10 

s stores the values T~ 1 (p, q) and the program for evaluating equation (3). A random access memory (RAM) 
11 is used for storing samples of the received signal R, intermediate results during the calculations, and the 
samples of the channel impulse response C. Central processor unit (CPU) 12 carries out the calculations 
under control of the program stored in ROM 10, and writes the results of the computation into the RAM 11. 
The CPU also controls the writing of the received signal samples into the RAM 11 and the reading of the 

ro channel impulse response values from the RAM 11. 

In the event that the system will be used for signals transmitted over different channels, a channel 
indication signal is also stored in the RAM 11. Different channel impulse response values for the different 
channels are computed and stored, and are read out according to the channel being used. 

The test signal T actually used in a computer implementation of the invention for video ghost 

75 cancellation was sine x = (sin x)/x. This function is shown in Fig. 3. The sine x function has the advantages 
of being well defined, containing adequate high frequency components to characterize the channel, and 
having sufficiently low values at its leading edge so that truncation errors are acceptable. Other test signals 
could be used as well. 

The matrix [T] is in Toeplitz form which has the property that the values of the elements on any 

20 particular diagonal are equal. In other words, a^ = a i + n , j + n - 

The computation of the inverse matrix [T] _1 can be simplified, and the amount of required memory in 
the receiver can be substantially reduced, by modifying it to truncate its size and place it in a cyclic 
Toeplitz form. How this is done is shown in Fig. 4A. 

In the example shown in Fig. 4A the sequence of {T(i)} has three members: {T(1), T(2), T(3)}. The test 

25 signal evolution matrix is a 6 x 6 matrix, and the received signal and channel vectors each have six 
elements. Assume in this example that T(2) has the largest value of the test signal samples. The matrix [T] 
is truncated so that the elements T(2) of each row fie on the main diagonal of the truncated matrix. The 
dashed lines in Fig. 4A indicate that the first row and the last column of the matrix [T] are deleted. This has 
the effect of eliminating the row and columns which are comprised of elements which are predominately 

30 zero. If the matrix [T] were to be inverted without deleting the row and column with predominately zero's 
and small valued elements, the resulting inverse matrix [T] _1 would have element values extending over a 
great magnitude range and make subsequent calculations involving the inverse matrix very difficult. 

The truncated matrix is put in cyclic form in the following manner. The element T(1) from the deleted 
row is placed in the first position of the last row, as shown by the arrow "a" in the drawing. Next, the 

35 element T(3) is "extended" diagonally to a position below the last column of the last row, as shown by the 
arrow "b", and then "reflected" to the first row as shown by arrow "c". The resulting matrix is Toeplitz, and 
it is also cyclic. Each column and each row of the matrix contains the test signal samples in cyclic order, 
i.e. if a column is traversed starting with the first test signal sample T(1), and if the traversal continues from 
the last element of the column to the first element of the column, then the test signal sample sequence will 

40 be traversed in order. The same holds true for each row. A cyclic Toeplitz matrix has the property that if it 
is invertible its inverse is also Toeplitz. In the rest of this description the matrix [T] will denote the test signal 
evolution matrix in cyclic Toeplitz form. 

The cyclic Toeplitz matrix derived from [T] as described above will in practice not always be invertible 
because the rank of the matrix [T] can be less than the dimension of the matrix. This occurs in television 

45 systems, for example, because the sampling frequency is substantially higher than the maximum video 
signal frequency, i.e. the video signal is oversampled. The NTSC picture signal has a color subcarrier of 
approximately 3.58 MHz which is sampled at approximately 14.3 MHz or four times the subcarrier 
frequency and twice the Nyquist sampling rate. This is done to allow the use of inexpensive filters while still 
avoiding aliasing of the sampled signal. 

so If the cyclic Toeplitz matrix [T] is subjected to an eigenvalue decomposition the resulting matrix [x] will 
be diagonal with the diagonal elements being the eigenvalues of the linear transformation represented by 
the matrix [T]. Some of the diagonal elements of [X] will be zero which is a consequence of the rank being 
less than the dimension of the matrix. 

If the matrix [T] is considered a linear transformation on some vector space, then the non-zero 

55 eigenvalues correspond to a sub-space (called the rank space), and the zero-valued eigenvalues correspond 
to another sub-space (called the null space). A property of these sub-spaces important to this invention is 
that they are disjoint; there is no orthogonal projection from one onto the other. Consequently, [T] will not 
be invertible. The mathematical details can be found in most texts on linear algebra, for example: Linear 
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to 



J5 



Algebra by Paul Shields, (Addison -Wesley 1964). 

In order to guarantee the invertibility of the matrix [T] and to limit the range of magnitudes of the 
inverse matrix elements, the applicant has discovered the following approximation technique. 

The test signal evolution matrix [T] is transformed by an eigenvalue decomposition transformation to 
obtain a matrix [X] of the eigenvalues of [T] as follows: 

[E] [T] [E]- 1 = [X] 

where [E] is the eigenvalue transformation and [X] is the matrix of eigenvalues 



[X] = 



25 
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30 



In general, the eigenvalues Xi to X, will be non-zero, and the eigenvalues X r+ i to X„ will be zero. The 
value of r, the number of non-zero eigenvalues, is called the rank of the matrix [T]. In order for a matrix to 
be invertible its rank must equal its dimension. This can be achieved by adding to [X] an augmenting matrix 
[A] which has non-zero values at the position of zero-valued eigenvalues of [X] such that 



35 



40 



45 



[X A ] = [X] + [A] = 



0 
X 2 



0 
0 

x r 



0 
0 



The inverse of the eigenvalue decomposition transformation [E] is then applied to [A] as follows 
so [W] = [E]"' [AJ [E]. 

This new matrix [W] has the property that every column of [W] is uncorrected with every column of [T]. An 
arbitrary one of the columns of [W] is then selected and added to the test signal sample sequence {T(i)} to 
form a new test signal sample sequence {T*(i)}. This new sequence {T*} is used to form a new test signal 
55 evolution matrix [T*] which is always invertible because it contains the non-zero eigenvalues of the 
augmentation matrix [A]. Additionally, when the matrix [T] is inverted and used to solve the equation 

C] = [TT' R] (4) 
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the uncorrelated components from W] will tend to average out. The channel impulse response sample 
sequence {C(i)} for the particular channel serves to determine the restoration filter coefficients. The input 
signal is filtered and the resulting output has channel-induced distortion diminished or removed completely. 

5 In order to compensate for variations in the channel impulse response which occur over time, the channel 
impulse response sample sequence {C(i)} can be repeatedly recomputed and the channel coefficients 
updated. Additionally, as mentioned in connection with the description of the circuit shown in Fig. 2, signals 
which have passed through different channels can be restored using the same filter by passing an 
appropriate set of channel impulse response samples to the filter. 

w Fig. 5 illustrates a general scheme for the filter configuration. A finite impulse response filter 20 receives 
the sequence {C(i)} of channel impulse response samples from either of the circuits shown in Figs. 1 and 2, 
and the filter coefficient values are set. The architecture of the filter 20 is arbitrary and the system designer 
may select among the numerous available designs based on the usual engineering trade-offs. 
The received signal sample sequence {X(i)} is applied to the input of a subtracter 22 as is the output of the 

75 filter 20. Here, the input signal X is any received signal and not just the test signal T. The output of the 
subtracter 22 is a sequence {Y(i)} of samples of the restored transmitted signal Y with the effects of the 
channel removed from the received signal X. 

As discussed above, the method and circuits according to the invention are completely general and rely 
on no prior assumptions regarding the channel characteristics, channel noise, the test signal and similar 

20 considerations. It exhibits high noise immunity and it is easy to implement. Accordingly, the particular 
circuits shown are exemplary, and the scope of the invention is defined by the following claims. 

Claims 

25 1. A system for communication channel identification, comprising: 

- a transmitter for transmitting a test signal T over a communication channel to be identified; 

- a receiver for receiving a signal R which is the test signal T after it has passed through the 
communication channel; 

characterized in that the system comprises calculating means for calculating a sequence of channel 
30 values 



Q 

{C(p)=52r- 1 (p,g)i?(g)} 

35 q=\ 

where T _1 (p, q) is representative for the (p, q)th element of the inverse of the test signal evolution 
matrix 

40 



45 



50 
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4. 
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6. 



{T(i)} is a sequence of values of the transmitted test signal and {R(i)} is a sequence of values of the 
received test signal. 

A system for communication channel identification according to Claim 1. characterized in that T~ 1 (p. q) 
is representative for the (p, q)th element of the inverse of a cyclic Toeplitz matrix derived from the text 
signal evolution matrix; an in that the system comprises means for discarding at least some of those 
elements form {C(m)} that have errors resulting from the use of the inverse matrix in Toeplitz form. 

A system for communication channel identification according to Claim 1 or 2, characterized in that T~ 1 - 
. (p,q) is representative for the (p,q)th element of the inverse of a matrix [V] derived from the test signal 
evolution matrix defined by a component being uncorrelated to each row of the signal evolution matrix 
added to the test signal evolution matrix. 

A system for communication channel identification according to Claim 3, characterized in that T~ 1 (p. q) 
is representative for the (p, q)th element of a matrix derived from matrix [T] by deleting a number of 
the upper rows of [V] and an equal number of the right most columns of [T*]. 

A system for communication channel identification according to on of the Claims 1 to 5, characterized 
in that the communication channel is a video channel, and the test signal T is a video test signal. 

A system for removing channel-induced distortion from signals, using a system according to one of the 
claims 1-5 characterized in that the system comprises means for filtering signals transmitted over the 
channel according to filter coefficients determined from the sequence {C(p)} of calculated channel 
values. 

Receiver suitable for use in a system as defined in one of the Claims 1 to 6. 



50 



A receiver according to Claim 7, characterized in that the receiver comprises 

- means for storing calculated channel impulse response samples for more that one channel; and 

- said calculating means is effective to calculate channel impulse response samples to permit 
multi-channel signal restoration. 



9. A method of communication channel identification, comprising: 

- transmitting a test signal T over a communication channel to be identified; 

55 - receiving a signal R which is the test signal T after it has passed through the communication 

channel; and 

- calculating a sequence of channel values 
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{c(p)=Y, T ~ { (P>3) R (<I)} 
4=i 



where T" 1 (p, q) is the (p, q)th element of the inverse of the test signal evolution matrix 
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{T(i)} is a sequence of values of the transmitted test signal and {R{i)} is a sequence of values of the 
30 received test signal. 

10. A method of communication channel identification according to Claim 1, further comprising: 

- writing the matrix [T] in cyclic Toeplitz form and inverting the matrix [T] to form an inverse matrix 
[T]~ 1 which is also in Toeplitz form; 

35 - calculating the sequence {C(p)} of channel values using the inverse matrix [T]" 1 in Toeplitz form; 

and 

- discarding at least some of those elements of {C(m)} that have errors resulting from the use of 
the inverse matrix [T]" 1 in Toeplitz form. 

40 



45 
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